Subwavelength microscope that uses frequency scanning for image reconstruction 
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A new principle of subwavelength imaging based on frequency scanning is considered. It is shown 
that it is possible to reconstruct the spatial profile of an external field exciting an array (or coupled 
arrays) of subwavelength-sized resonant particles with a frequency scan over the whole band of 
resonating array modes. During the scan it is enough to measure and store the values of the near field 
at one or at most two points. After the scan the distribution of the near field can be reconstructed 
with simple post-processing. The proposed near-field microscope has no moving parts. 



I. INTRODUCTION 

The maximum achievable resolution of any conven- 
tional lens, the function of which is based on focusing of 
electromagnetic radiation, is the wavelength of the used 
radiation (the limit that is also known as the diffrac- 
tion limit or the Rayleigh limit). The resolution limit 
is based on the fact that the evanescent waves, which 
carry the information about the subwavelength proper- 
ties of the source, are lost due to their exponential atten- 
uation, and only the propagating waves are restored in 
the image plane by a focusing lens. The resolution of an 
imaging device can be greatly improved if this device is 
able to react to the evanescent part of the spatial spec- 
trum of the incident field. There are several known tech- 
niques which allow subwavelength-resolution imaging: 
superlenses based on the use of materials with negative 
parameters, 1,2 superlenses based on phase-conjugation 
(non-linear three- wave mixing) , 3 arrays of small resonant 
particles. 7-11 However, the only mean to actually mea- 
sure the image-plane near-field distributions at subwave- 
length scale is the scanning near-field microscope which 
uses a small moving probe. Recently, hyperlenses which 
transform evanescent modes into propagating ones have 
been proposed 4-6 with the goal to develop subwavelength 
microscopes. These devices allow image detection by a 
stationary system like a charge-coupled device (CCD) 
matrix, however, the manufacturing technology for hy- 
per lens structures needs further development. 

In this paper we introduce a near-field subwavelength 
imaging device that is based on frequency scanning with 
following post-processing of measured data. The field is 
measured only at one or two points in space, and the role 
of sensor is played by an electrically dense grid of small 
resonant particles. The device does not use any moving 
probes, allowing for very fast subwavelength imaging. In 
this device the conventional measuring of near fields at 
many points in space is replaced by measuring resonantly 
enhanced near fields at many frequency points, which 
allows us to calculate the spatial field distribution. The 
maximum spatial resolution of the device is limited by 



the period of the sensing grid and depends on the quality 
factor of grid particles as well on the number of particles. 

The present method is a development of the approach 
of our previous paper 7 which introduced a system of two 
coupled resonant grids or arrays of small inclusions as 
a superlens capable for resonant amplification of evanes- 
cent fields and creation of images with subwavelength 
detail. That idea was further extended for enlarging 
superlenses. 8 Possible realization in the optical region 
was also demonstrated. 9,10 Experimental confirmations 
in the microwave region were made using grids of small 
resonant electric dipoles 7,8 and magnetic dipoles (split 
rings). 11 Grids of resonant particles get strongly ex- 
cited when the transversal wavenumber k t of an incident 
evanescent wave matches the propagation factor of grid's 
surface wave. We showed both theoretically and experi- 
mentally how this effect can be used for resonant ampli- 
fication of evanescent fields. Because of the resonance, 
the amplitude of the evanescent wave in the image plane 
of a dual-grid structure 7 becomes equal to the amplitude 
of the evanescent incident wave in the source plane. 

Matching of the propagation factors takes place at a 
certain frequency. The grid of densely packed resonant 
particles can be modeled by its grid impedance Z g (ui,kt), 
which connects the surface-averaged tangential electric 
and magnetic fields. 12 The grid impedance depends on 
the frequency uj and the propagation factor along the 
grid plane fc t . The dispersion equation for the surface 
waves on an impedance grid in free space can be written 
as 

Z g {w,kt) + =0, (1) 

where Zq(lj, k t ) is the wave impedance of the correspond- 
ing free-space plane wave: 

Jh.2 _ ju2 

Z («, fct) = -> TM modes, (2) 

fc 

Zo(u>, k t ) = T7 =, TE modes, (3) 

y/ K — k t 

where 770 = \/ 77 an d fco = <jJ^hq£q. The solutions of 



2 



this equation are pairs (u>, k t ) that lie on grid's dispersion 
curve. If the grid is periodic with the period d -C Ao (Ao 
is the free-space wavelength), its dispersion curve usually 
looks like it is shown in Fig. 1(a). 
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FIG. 1: a) An example dispersion curve of a periodic struc- 
ture. On the horizontal axis q — k t d. b) Transmission trough 
a pair of resonating arrays at a given frequency as a function 
of q. 

We see that for an operating frequency lu G 
[w m i n , w max ] there is only a single resonating value for the 
propagation constant k t . Scanning the frequency from 
w m i„ up to w m ai we can sequentially excite the surface 
modes of all possible kt ranging from to fc fmax = n/d. 
The amplitude of an excited mode is proportional to the 
amplitude of the corresponding spectral component of 
the incident field. Knowing the spatial profile of every 
surface mode (from the theory or from an initial calibra- 
tion measurement) we can find the amplitude of the mode 
from a field measurement done only at a single point in 
the image plane (or at most at a couple of discrete points, 
to account for degenerate cases of zero field of certain 
modes at the measurement point). In a finite-size grid 
the spectrum of surface modes is discrete. In this case 
the field measurement is done at discrete frequencies. 

It is worth noting that for open structures (in con- 
trast to setups confined inside a closed impenetrable cav- 
ity) the surface modes with < k t < fco <C ir/d are 
leaky modes, i.e., they radiate into surrounding space. 
However, the operation of the microscope is based on 
the use of eigenmodes with high propagation constants 
(kt 3> fco), which are strongly bound to the grid sur- 
face and radiate only at inhomogeneities and at the grid 
ends. Moreover, in some resonant systems such as the 
metasolenoid 13 even the modes with small values of k t 
produce a strong resonant response because of a high Q- 
factor of its resonances (Q is of the order of 10 3 ). 



In a system of two coupled resonant grids 7 the am- 
plitude of the field in the image plane can be made the 
same as in the source plane so that the system works as 
a "perfect lens" 1 at every resonant frequency of the sys- 
tem. This allows us to find the spatial field distribution 
of the unknown source by weighted integration of mea- 
sured values, as is explained in Section II. In principle 
one can use just a single resonating grid. The only diffi- 
culty here is that the relation between the amplitude of 
the excited surface mode and the unknown external field 
is more complex. The single-grid microscope is explained 
in Section III. 



II. PRINCIPLE OF SPATIAL IMAGING BY 
FREQUENCY SCANNING 

For clarity, let us consider a lens composed of a couple 
of arrays of small resonating electric dipoles (oriented 
along the a;-axis) placed in between two metal screens: 
the same system as we used in our measurements. 7 The 
dispersion curves of the arrays look like in Fig. 1(a) and 
the field is measured at the the middle point of the im- 
age plane at x = 0. To cope with the situation when the 
x-component of the electric field vanishes at this point 
for a certain mode it is necessary to have a couple of sen- 
sors there: one measures E x and the other measures the 
orthogonal component E z (the longitudinal component). 
We will also assume that the arrays consist of many par- 
ticles, so that the spectrum of its modes is practically 
continuous. An example geometry of the structure is 
shown in Fig. 2. 
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FIG. 2: The side view of a possible lens structure. Two ar- 
rays of small particles are positioned along the x-axis. The 
screening box is open at one of its sides (at the source plane) . 

Let us assume that the object under study is located 
in the source plane and it can be excited at the frequen- 
cies in the range [u; m i n , w max ] and that the properties of 
the source itself do not change much when we scan the 
frequency in this range. Practically this means that we 
have to provide as narrow frequency band of the disper- 
sion curve as possible. In terms of the design of the grids 
this leads to the known conditions: 7 the particles must be 
weakly interacting and be high-Q resonators themselves. 
At microwave frequencies, electrically small copper- wire 
meanders are good enough with enough narrow resonant 
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band. 7 In the optical range plasmonic metal spheres can 
be used. 9 ' 10 

If the source is passive (for example, it is a collection of 
small pieces of metal) it can be excited with an antenna 
that illuminates it with a plane wave of the necessary fre- 
quency. Because the structures we deal with are most ef- 
fectively excited with evanescent waves, the illuminating 
source will not influence the operation of the device (or 
at least its influence is predictable and can be taken into 
account in the post-processing stage). It can be a good 
engineering task to figure out the optimal disposition of 
the illuminating source and the arrays. In the particu- 
lar example considered here (Fig. 2), it is assumed that 
the source electric field is polarized along x and traveling 
modes are not supported inside the screening box. 



For a given frequency lu G [lu 



mm; ^max 



the transmis- 
sion coefficient through the lens as a function of k x will 
look like in Fig. 1(b). The characteristic width of the 
spike is physically determined by the characteristic size 
L (the array length) of the lens: Ak x rs tt/L. In principle 
the spike can be made as narrow as required. The field 
sensor will measure a kind of average amplitude of the 
modes with k XrcB — Ak x /2 < \k x \ < k Xrcs + Ak x /2. 

The device operates as follows. We scan the frequen- 
cies from w m i n up to w max and store the measured elec- 
tric field complex amplitudes at every frequency. In the 
post-processing we use the known spatial profiles of the 
modes to restore the actual field distribution from the 
measured field values at the source plane. This is possi- 
ble via simple integration, because the amplitudes of the 
fields of every mode are the same at the image and source 
planes (neglecting dissipation in the grids). 7 Denoting 
for even modes: A cvcn = E^ mhc , and for odd modes: 
A odd = -(k z (u)/k x (uj))EP mhc , the field distribution in 
the source plane can be expressed as a sum of all modes: 



E x (x) = - 

IT 



A cvcn {uj)cos(k x (Lu)x) 



+ A oM {u})sm(k x (u)x) 



dk x (lo) 



du. (4) 



The coefficient L/ir accounts for the spectral density of 
modes. 

As it was said above, in a finite-size grid the spectrum 
of surface modes is discrete. In this case the field mea- 
surement is better to be done at discrete frequencies. We 
will consider such a system in the next section. 




FIG. 3: A linear array of plasmonic nanoparticles as a sub- 
wavelength microscope sensor. The line shows the electric- 
field profile to be determined by the microscope. The probe 
is stationary. 



the array and that the array period is d. The electric 
dipole moments pk of the array particles can be found 
from the following system of linear equations: 



N 



_1 (w)pfc 



/3 n -k(co)p n = E\ 

71=1, n^k 



ext 
k ' 



(5) 



where E^, xt is the external electric field at the position of 
the particle number k, a{ui) is the dipole polarizability, 
and P n -k(u>) accounts for interaction in a pair of dipoles 
with indices n and k. Our goal is to find the spatial dis- 
tribution of the external field (that is, find vector E^ xt ) 
from measured field at one point but at several frequen- 
cies. 

For dipoles orthogonal to the line of the array 
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For dipoles parallel to the line of the array 
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The polarizability of a small resonant dipole can be ap- 
proximated as 



1 



a(uj) 



1 



LU Q 



6irc 2 ' 



(8) 



where ao is the static polarizability, c^o is the resonant 
frequency, Q is the quality factor of the particle and the 
last term accounts for the radiation loss. 12 

The system of equations (5) can be written in matrix 
form 



III. SINGLE ARRAY OF RESONANT 
PARTICLES AS A NEAR-FIELD SENSOR 

Let us consider a linear array of resonating particles, 
Fig. 3. The particles are electrically small passive in- 
clusions that we model by electric dipoles with known 
polarizabilities. We assume that there are N particles in 



A(u) -p{lu) = E° 



(9) 



where A(uS) is the system matrix, the elements of which 
can be calculated from the above formulas; p(uS) and E 
are vectors of length N . The eigenmodes of the array can 
be found from the equation 



detA(w) = 0. 



(10) 
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There are N modes in total. In a system with non-zero 
loss the modal frequencies defined in this manner are 
complex quantities with the imaginary part representing 
the decay of oscillations. In our problem, however, the 
frequency is defined by the applied external field and is 
real. Therefore, the array resonances are found as the 
solutions of equation 



Rc{det A(oj k )} = 0, 



(11) 



where all Lo k are real. Note that the matrix A{ui) is almost 
singular at the frequencies defined by (11). Therefore, its 
inverse can be approximated as 



A ( Wfc ) « \-\u k )V(w k ), 



(12) 



where X(u>h) is the eigenvalue of A{uu) with the smallest 

magnitude at ui — > ui k and V(u)k) is the projection opera- 
tor for the corresponding eigenvector. Hence, the dipolc 
moments at the resonances are 



(13) 



From here we see that at the resonances the array acts 
effectively as a filter extracting a single spatial harmonic 
from the external field. 

At the resonances the induced dipolc moments arc 
maximal and so is the electric field at the probes that 
we place near the array. Therefore, it makes sense to 
measure these fields at resonances. At these frequencies 
u>k the field at the location of a probe can be found as 



N 



(14) 



n=l 



where c n (co k ) are the coefficients describing coupling of 
the probe to each dipolc in the array. 
We can rewrite (14) in matrix form 



E(uj k ) = c T (uj k ) ■ p(uj k ) = c T (uj k ) ■ A (u> k ) ■ if 
or, for all N frequencies at once 



\c T (uj n ) ■ A (o> n ), 



(15) 



(16) 



where the rows in the central matrix are products of c T 

and A at all the resonant frequencies. Let us recall 
that we assume that the distribution of the external field 
does not change when we scan the frequencies in the 
range [uji,u) N ]. 

Finally we can solve (16) for the distribution of the 
external field: 



'Ef 



TT'CXt 

1 N 



\c T {uj n )-A (w„), 




(17) 



For this solution to exist the matrix in (16) must be non- 
singular. This can be achieved either by placing a probe 
at a non-symmetric location with respect to the dipole 
array or by using more than one probe and intermixing 
their responses in (16). 



IV. NUMERICAL EXAMPLE 

As an example let us calculate the response of a lin- 
ear array of orthogonal dipolcs and then restore the dis- 
tribution of the external field using the method intro- 
duced above. We choose the following parameters for 
the dipole array: the period is such that kod = 0.1 at 
the dipole resonance, the number of dipoles in the array 
N = 21, the quality factor of the particle itself (with- 
out radiation loss which is accounted for separately) is 
Q = 10 3 , the static dipole polarizability is such that 
a /(eod 3 ) = 0.5. This value of the polarizability is rather 
high. Fo r a metal sphere it corresponds to the radius of 
a = d-^/0.5/(47r) rts 0.34 d. Note, however, that reso- 
nant particles of complex shape or plasmonic particles 
can provide high values of polarizability even if their size 
is electrically quite small. 

In Fig. 4 the resonant frequencies of the array are plot- 
ted as function of their index varying from 1 to 21. These 
frequencies are found from (11). 
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FIG. 4: Resonant frequencies of the dipole array. 

The probe is located close to the array at the distance 
of half a period from the array line. For the probe to 
be able to react to both symmetric and antisymmetric 
modes the probe must not be placed in a symmetric po- 
sition with respect to the array. To find the optimal 
position we investigate how the condition number of the 
matrix (16) changes when we move the probe along the 
array. This is shown in Fig. 5. 

The optimal position is found to be in front of the 
second dipole if counted from the middle of the array. 
At this point the condition number is about 10. This 
means that the measurement errors or any other noise 
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FIG. 5: Condition number of the system matrix as a function 
of the probe location. Coordinate x is along the array with 
the origin at the middle dipole. 
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FIG. 7: The profiles of the external field (solid lines) and the 
restored field (symbols). 



of relative amplitude of 0.1 will destroy the restoration 
process. For practical purposes the measurement error 
must not exceed 10 _3 -10~ 2 . 

The higher is Q the lower is the condition number. The 
same holds for the value of a . From the other hand, the 
condition number grows with increasing the number of 
dipolcs in the array. For an array of about 50 particles 
to obtain the condition number of 20 one has to provide 
Q ~ 10000 which is hardly achievable. Therefore, with 
large arrays it would be more practical to employ several 
probes placed at carefully selected positions. 



in Fig. 6. The probe measures the field only at the ar- 
ray resonant frequencies. In Fig. 7 the restored spatial 
profile of the external field is compared with the origi- 
nal one given by (18). To show the effect of noise we 
have added a normally distributed noise of the relative 
variance of 10~ 2 to the measured probe field. If there 
is no noise, the restoration is exact, that is, the values 
of the incident field at the positions of each particle in 
the array are determined exactly. This means that the 
microscope resolves spatial details down to the scale of 
one array period. 
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FIG. 6: Electric field amplitude at the location of the probe 
at different resonant frequencies. 

In the following field restoration example we assume 
the following distribution of the external field: 



1 



1 



l + (k-9f J l + (fc-ll) 2 ' 



(18) 



With such an external field in action the amplitude of 
the electric field at the location of the probe is shown 



V. CONCLUSIONS 

In this paper we have presented a near-field imag- 
ing principle which allows us to determine spatial 
distributions of electric or magnetic field with sub- 
wavelength resolution using observations at several fre- 
quencies within a narrow frequency range of sensor's res- 
onances. The field is measured only at one or a few se- 
lected points behind the sensing array of resonant parti- 
cles. The limit for the image resolution is determined 
by the grid period. We have shown that for reliable 
measurements it is desirable to use resonant particles 
with high quality factors and high polarizability val- 
ues. The evanescent fields emanating from the object 
under study are effectively enhanced due to resonant ex- 
citation of the sensor particles. The necessary range 
of the frequency scan is reduced if the array particles 
are weakly coupled to each other. Since the device has 
no moving parts and there are known means to change 
the frequency of the probing electromagnetic wave very 
quickly, this method offers a possibility to make obser- 
vations of complex sub- wavelength objects very fast and 
in a non-invasive way. For microwave applications, the 
sensor can be formed by electrically small metal par- 
ticles, e.g., in form of a meander 7 or as dense pack- 
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ages of split rings (metasolenoids). 13 For optical appli- 
cations, the sensor can be potentially realized as grids 
of plasmonic nanoparticles. 9 ' 10 Another possibility is to 
use grids of resonant nanocavities in thin metal layers. 14 
The last case has the advantage of reduced coupling be- 



tween resonant voids due to field decay in the metal layer. 
We are now studying dispersion properties of arrays of 
nanovoids in metal near an interface with free space with 
some promising results. 
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